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ABSTRACT 

We study the initiation of the detonation in the gravitationally confined det- 
onation (GCD) model of Type la supernovae (SNe la). In this model, ignition 
occurs at one or several off-center points, resulting in a burning bubble of hot 
ash that rises rapidly, breaks through the surface of the star, and collides at 
a point on the stellar surface opposite the breakout, producing a high- velocity 
inwardly directed flow. Initiation of the detonation occurs spontaneously in a 
region where the length scale of the temperature gradient extending from the 
flow (in which carbon burning is already occurring) into unburned fuel is com- 
mensurate to the range of critical length scales which have been derived from ID 
simulations that resolve the initiation of a detonation. By increasing the max- 
imum resolution in a truncated cone that encompasses this region, beginning 
somewhat before initiation of the detonation occurs, we successfully simulate in 
situ the first gradient-initiated detonation in a whole-star simulation. The deto- 
nation emerges when a compression wave overruns a pocket of fuel situated in a 
Kelvin- Helmholtz cusp at the leading edge of the inwardly directed jet of burn- 
ing carbon. The compression wave pre-conditions the temperature in the fuel in 
such a way that the Zel'dovich gradient mechanism can operate and a detona- 
tion ensues. We explore the dependence of the length scale of the temperature 
gradient on spatial resolution and discuss the implications for the robustness of 
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this detonation mechanism. We find that the time and the location at which 
initiation of the detonation occurs varies with resolution. In particular, initia- 
tion of a detonation had not yet occurred in our highest resolution simulation 
by the time we ended the simulation because of the computational demand it 
required. However, it may detonate later. We suggest that the turbulent shear 
layer surrounding the inwardly directed jet provides the most favorable physical 
conditions, and therefore the most likely location, for initiation of a detonation 
in the GCD model. 

Subject headings: hydrodynamics — shock waves — supernovae: general — white 
dwarfs 



Introduction 



Type la supernovae (SNe la) are among the most luminous optical transients. They con- 



tribu te significantly to the non-primordial elements in the universe, especially iron (ITimmes et al. 



19951 ). A large fraction of SNe la exhibit a correlation between peak luminosity and rate 
of decline that can be used to make them standard candles, and therefore excellent dis- 



tance indicato r s and probes of cosmology [see, e.g. the reviews by iLeibundgutl (120011 ) and 
Frieman et al.l (120001 )]. While S Ne la are thought to b e the result of the thermonuclear incin- 
eration of white dwarfs (WDs) (IHoyle &: Fowlerlll960l). the explosi on m echanism is not fully 



under stood [see e.g., the reviews by iHillebrandt &: Niemeyerl (120001 ) and lPodsiadlowski et al. 

( boosh ]. 



Most astronomers favor a scenario in which the exploding star is a single, near Chandrasekhar- 
mass WD. Consequently, multi-dimensional full-star simulations of the explosion of near 
Chandrasekhar-mass CO white dwarfs have been performed by several research groups. 
Thes e models posit the ignition of a turbulent deflagration flame in the c ore of the WD 



Garcia-Senz fc Wooslevlll995l : iHoflich fc Steinll2002l ; Kuhlen et alJl2006h . which eventu 



(e.g. 

ally gives rise to a supersonic detonation wave. 

In the deflagration to detonation transition (DDT) model, the detonation is postulated 
to occur when the deflagration flame reaches densities low enough that the Gibson length l@, 
defined as the scale where the turbulent velocities equal the laminar flame spe ed, becomes 
smaller than the flame width 5 (IKhokhlov et al.lll997l ; lNiemeyer &: Woosleylll997l ). When this 
criterion is met, turbulence tears the flame apart, producing distributed burning. Simulations 
of SNe la that explode via DDT have been extensively studied in multiple dimensions (e.g . 
Arnett fc Livn3l 1 994al fbL lLivnelll999l : bamezo et albooi 120051 : iGolombek &: Niemeverll2005l : 
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Ropke fc Niemeyerl 120071 ; iBravo &: G arcfa-Sena 120081 ) . 



The pulsational delayed detonation model (PDD) posits a centrally ignited deflagra- 



Ivanova et al. 


1974; 


Khokhlov 


1991) 



onation (PRD) model, in which the initial deflagration releases energy and expands the 
star but quenches before it becomes gravitationally unbound. The detonation in these 
models is postulated to form from an accretion shock during the re-collapse of the star 
after the deflagration ash has risen and fresh fuel occupies again the center of the WD 



(jDunina-Barkovskaya et al.ll200ll ; IBravo &; Garcfa-Senzll2005l . 120061 ) . 



In the gravitationally confined detonation model (GCD), the deflagration flame ignites 



at one or more points s l ightly offset from the center of the WD (IPlewa et al 



20071 : iRopke et all 120071 : iTownslev et all 120071 : IJordan et all 120081 : Meakin et al 



2004 



Plewa 



2009). As 



the bubble grows and rises, the Rayleigh- Taylor instability causes rapid growth of the flame 



Khokhlov 


1995: 


Zhang et al. 


2007) 



erates, burning on the order of a fe w percent of the m ass of the star by the time it breaks 



through the surface of the WD (e.g. iLivne et al.l 120051 ) . Owing to the large surface gravity 



of the WD, the bulk of the ash in the bubble remains gravitationally bound. The cold ash 
flows over the surface of the star and converge at the opposite point on the stellar surface 
from where breakout occurred. The converging flow of cold ash pushes unburned fuel in the 
surface layers of the WD ahead of it, compressing it, and increasing its temperature and 
pressure enough that the carbon in it begins to burn. Inwardly and outwardly directed jets 
form at the s tagnation poi n t, dri ving the inwardly directed jet to high densities [ see the 
discussion in iMeakin et al.l (120091 )]. The inwardly directed jet is unstable to shear insta- 
bilities and Kelvin-Helmholtz rolls are observed to form. Additionally, strong sound waves 
originating from the surface flow traverse the jet. 

It is in this environment that the detonation in the GCD model is thought to occur. To 
determine whether initiation of a detonation occurs under these conditions is, however not 
straightforward, primarily because the length scales pertinent to the detonation initiation 
are typically much smaller than the highest resolution possible in 2D, and even more so in 
3D, whole-star simulations. In this context, past simulations of the GCD scenario can be 
grouped into two categories: 



1. Those that limit nuclear burning to the flame model and don't simulate the detonation 
phase, including in itiation of the detonation (IRopke et al.l 120071 ; ITownslev et al.l 120071 : 
Jordan et~al1l2008h . 
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2. Those in whi ch the detonation emerges in situ through nuclea r burning outside the 
flame model Jpiewa et alibooi IPlewal I2OO7I : iMeakin et aDl2009h . 



Simulations in the first category only infer the successful initiation (or failure) of the 
detonation based on a comparison of the temperatures and densities reached in the full-star 
simulation with those for t he "critical detonation condi tions" derived from ID reactive hydro- 
dynamics simulations (e.g. iNiemeyer fc Woosleyl 119971 ) at higher resolution. This procedure 
is inherently uncertain, due to the sensitivity of the detonation con ditions to the functiona l 
form of the unresolved and therefore unknown temperature profile (jSeitenzahl et al.ll2009l ). 
The practice of using the temperatures and corresponding densities in a simulation without 
distributed burning (i.e. in which nuclear burning is not allowed to occur outside the flame 
model) suffers from additional limitations. First, it ignores the feedback of nuclear burning 
on the hydrodynamics. The energy liberated during carbon burning in the jet alters the 
flow and may produce a temperature structure more conducive for initiation of a detonation. 
Second, not keeping track of nuclear reactions takes away the ability to decide whether the 
temperatures and densities obtained in the purely hydrodynamical flow in fact occur in det- 
onatable fuel, or the fuel has already been depleted due to the nuclear burning that occurs 
during the process of compressional heating. 

It is our belief that direct initiation of a detonat ion is hard to achieve in a CO WD 
environment (INomoto et al.lll976l ; iMazurek et al.lll977) and that spontaneou s initiation via 



the gradient (SWACER) m e chanism (jZel'dovich et al 



19701 ; iLee et al.l Il978l ) is more likely 



telinnikov fc Khoklovl [l986l : iBlinnikov fc Khokhlovill987h . GCD simulations in the second 



category allow for the detonation to form explicitly in situ, and have paid some attention to 
the details of the initiation mechanism. However, no GCD simulations to date have spatially 
resolved the temperature gradients that are key to initiation of the detonation. 

In section [2] we describe the setup for the simulations reported in this paper, and the 
method we use to resolve the temperature gradients in the region where the detonation is 
expected to occur and so shed more light onto the details of the detonation mechanism in 
GCD. In section [3] we present our results for a set of four simulations with different maximum 
resolutions. Section H] compares our results with earlier work and discusses the implications 
of our results. Section [5] summarizes our conclusions. 



2. Simulations 



In this paper we present 2D axisymmetric GCD simulations perfo rmed with the flas h 
Code. The way in which we carry out the simulations is identical to IMeakin et al.l (120091 ). 
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who extend the ITownslev et al.l (120071) de flagration models t o the d etonation phase and use 
the energetics of ICalder et al.l (120071 ) and ISeitenzahl et al.l (j2009bl ). The initial model is a 
cold (T = 4 x 10 7 K), isothermal WD with a mass of 1.365 M Q , which has a central density 
of p c = 2.2 x 10 9 g cm -3 . The composition throughout the WD is equal parts 12 C and 16 
by mass. The simulations start with a single ignition bubble with radius 16 km offset from 
the center of the star by 40 km offset. In this single ignition point scenario, the center of the 
bubble and the center of the star define a symmetry axis. The bubble rises as it grows and 
breaks through the surface of the star at what we will refer to as the "North Pole." Most of 
the ash in the bubble remains gravitationally bound to the white dwarf. Now rarefied and 
cold, the ash spreads out from the North Pole and flows over the stellar surface, compressing 
unburned fuel near the the opposite point on the stellar surface (which we will refer to as 
the "South Pole"). A conical accretion shock forms 300-400 km away from the symmetry 
axis (see Fig. [[]). The pressure in the compressed region exceeds the hydrostatic pressure in 
both radial directions, producing outwardly and inwardly directed jets (see Fig. [2]). 



2.1. Jet morphology 



The hydrodynamical flow in the collision region is complex: Compression waves emanate 
from the head of the inwardly moving jet into the star; the (locally) supersonic flow in the 
jet gives rise to internal shocks; strong sound waves from the i nteract i on of the surface flow 
with the underlying star traverse the South Pole region (see iPlewal (120071 ) and Fig. [3] ). 
Furthermore, the boundary between the jet and the surrounding material is unstable to the 
generation of turbulence via the Kelvin-Helmholtz (K-H) instability and "rolls" characteristic 
of the K-H instability form and broaden the interface. It is in this dynamic, hot, and 
turbulent environment that conditions may occur that lead to the initiation of a detonation. 

At the accretion shock, the press ure of the shocked m aterial is balanced by the ram- 
pressure of the inflowing material (c.f. Meakin et al.l I2OO9I ) . Compression of the fuel heats 
it to temperatures high enough for carbon burning to occur. This raises the pressure even 
further above the hydrostatic value. Save for a narrow channel o f cold material, carbon, as 
measured by the progress variable 0i (c.f. ITownslev et al.l 120071 ) . has been depleted in the 
inwardly movin g jet (see fig. IHleft). O n the other hand, oxygen (as measured by the progress 



variable 



(c.f. 



Townsley et al.l 120071 ) ) . is only locally depleted (mostly near the symmetry 



axis) and remains largely unburned in the jet (see Fig. [H right panel). The hot material in 
the jet is driven inward up the density gradient while continuing to consuming carbon. 
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2.2. Higher resolution truncated cone 

A goal of this paper is to determine whether the conditions near the head of the in- 
wardly directed jet in the GCD model are conducive to the initiation of a detonation via 
the Zel'dovich gradient mechanism or another mechanism (e.g., SDT; see the discussion in 
section H]). One-dimensional simulations show that, while the critical length scale for sponta- 
neous initiation of a detonation depend s on the temperature p rofile, at a density of 1.0 x 10 7 



g cm 3 it is on the order of a few km (jSeitenzahl et al.ll2009l ). Therefore, simulations with 



a resolution at least this fine are necessary in order to determine whether initiation of a 
detonation occurs via the gradient mechanism 

In spite of performing our calculations only in 2D and making use of the adaptive mesh 
refinement (AMR) capabilities of the flash code, full-star simulations at sub-km resolution 
were prohibitively expensive given the computational resources available to us. We therefore 
decided on a compromise, and introduced a higher resolution truncated cone covering the 
region of the jet where the detonation is expected to occur. By spatially constraining the 
highest resolution volume in this way, we limited the total number of computational cells 
while retaining the capability to resolve the length scales on which critical temperature 
gradients were expected to occur. 

The deflagration phase of all of the simulations presented here was identical and per- 
formed at 4 km resolution up to t — 2.2124 s. At this time a higher resolution truncated 
cone with (half) opening angle a ~ 10° was imposed on a sub-domain extending in radius 
from 1600 km to 2600 km (see fig. [5]). We chose the parameters of the highly resolved 
region so that it covered the compression region ahead of the inwardly moving jet at the 
South Pole where we expected initiation of the detonation to occur, assuring that this region 
was maximally refined. We carried out four simulations: one in which the resolution in the 
collision region ahead of the inwardly moving jet remained 4 km, and three others in which 
the maximum resolution in this region was 500 m, 250 m and 125 m. 



3. Results 

All but the highest resolution case (125 m) had detonated by the physical time at which 
we stopped the simulations. Table [1] gives the physical conditions (i.e., the time, location, 
peak and background temperatures, and fuel density) in the environment where initiation of 
the detonation wave occurs in the different cases. In the 500 m and 250 m cases, initiation 
of the detonation is triggered by the interaction of a hydrodynamical compression wave with 
the leading edge of the jet. This makes the definition of fuel density somewhat ambiguous, 
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as one could refer to the state on either side of the compression front. Tabled] therefore lists 
both the uncompressed fuel density (left) and the compressed fuel density (right) in these 
two cases. 

All of the detonations that occur in this study are due to the gradient mechanism; i.e., 
they emerge from a gradient in induction times. We calculate the isochoric induction time, 
T i,cv(P8, T 9 , Xq), defined here as the amount of time that passes until 90% of the initial mass 
fraction of 12 C, Xq, is consumed in a constant volume self heating network calculation, using 
the formula, 



r ijCU = 1.18 x 1(T 5 s (X^ps)- 1 ' 85 ^ - 0.206)" 7 ' 7 [1 + 1178(T 9 - 0.206) 



-1.85 



7.7 r 



-7.7l 



X 



This expression is a fit made by iDursi fc Timmed (120061 ) to a large number of 13-nuclide a- 
chain nuclear reaction network calculations, which the authors found was accurate to better 
than a factor of 5 for the densities and temperatures of interest here. For an in depth 
discussion of the spontaneous initiatio n of detonations via the gradient mechanism and the 
critical length scale for detonation, see ISeitenzahl et al.l ( 120091 ). 



Since the mechanism(s) leading to initiation of a detonation vary in the four cases we 
carried out, we discuss them individually. 



3.1. 4.0 km resolution case 

In this case (which is the lowest resolution simulation we carried out), the detonation 
wave occurs on the symmetry axis (see Fig. [6]) when the head of the jet reaches a fuel 
density of p ~ 1.1 x 10 7 g cm -3 . It occurs ~ 0.06 s earlier than in the 500 m and 250 m 
cases discussed below (see Tabled]). The detonation wave spreads radially outward from the 
initiation site. Ignition of oxygen and the products of carbon burning occurs concurrently 
with initiation of the the detonation wave (see fig. [6] right). 

Initiation of the detonation occurs via the gradient mechanism as a result of the tem- 
perature gradient that extends into unburned fuel ahead of the jet, which consists of carbon 
burning ash and is hot. The length scale of the temperature gradient is ~ 3 — 4 compu- 
tational zones, the number of zones over which flash spreads the temperature jump at 
the head of the jet, and thus remains constant. Since the critical length scale for initiation 
of a detonation decreases as the density increases, the detonation occurs when the density 
becomes high enough that the length scale of the induction time gradient (which for the 
approximately constant density conditions that exist locally in the unburned fuel is in one- 
to-one correspondence with the temperature gradient), becomes smaller than the critical 
length scale for initiation of a detonation. We consider the initiation of the detonation in 
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this case to be an artifact of the low resolution and numerical diffusion, which lead to a 
much shallower temperature gradient, and hence larger pre-heated region, than in the higher 
resolution cases. 



3.2. 500 m resolution case 

In this case, the detonation occurs off axis (see Table [1] and Figs. H]and[7j). It occurs 
~ 0.06 s later than in the 4 km case and at a density of ~ p — 1.2 x I0 7 g cm -3 , which 
is higher than in the 4 km case. The reason is that the length scale of the temperature 
gradients at the leading edge of the jet are still ~ 3 — 4 computational zones; i.e., the are 
an eighth of what they were in the 4 km case. Consequently, the temperature gradient at 
the time and location where the detonation was launched in the 4 km case are now too 
steep. However, t he cr itical temperature gradient is steeper at high densities (see Table 4 in 



Seitenzahl et al.l (120091 )). Detonation at the higher resolution is merely delayed until the jet 



has reached densities high enough that the critical length scale is satisfied. 

The detonation emerges behind a compression front that moves more or less transverse 
to the tip of the jet (see Fig. [7] left). The compression wave is not strong enough, i.e., the 
temperature behind the wave is not high enough (see Fig. [7J right), for direct initiation of 
a detonation. The wave, however, plays a crucial role in the initiation. It compresses and 
heats the carbon in the transition layer at the head of the jet, leaving behind an induction 
time gradient which eventually becomes shallow enough to lead to a detonation as the jet 
moves to higher density (see Fig. [8]). 



3.3. 250 m resolution case 



In this case, the detonation forms off axis when the leading edge of the jet has advanced 
to a density of ~ p = 1.2 x 10 7 g cm -3 (see Table [1] and Fig. [9]) in a location close to where 
the 500 m case detonated. The detonation occurs only fh 0.002 s earlier than in the 500 
m case. While both the 500 m and the 250 m simulations rely on the moving compression 
front for initiation, and detonate in a similar location, the mechanism is slightly different in 
the two cases. The 250 m model detonates when the temperature gradient at the head of 
the jet is broadened over a few zones extending into unburned fuel (see Fig. middle left). 
The shallow temperature gradient is established when a compression front and an associated 
reflected wave interact with a Kelvin-Helmholtz cusp (see E]), which is not present in the 
lower resolution simulations. The high temperature extends into a region of unprocessed 
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fuel and a relatively shallow gradient in the isochoric (constant volume) induction time is set 
up (see Fig. [8] middle right). The coherent energy release due to carbon burning leads to the 
build up of a shock with high enough pressure (see fig. [TU]) that in its wake oxygen burning 
commences along a line along the original induction time gradient (see Fig. fTT]) . Once oxygen 
has ignited, the detonation begins to propagate radially outward from the initiation site. 



3.4. 125 m resolution case 

This case initially resembles the 250 m case. The same hydrodynamical wave that 
triggered detonation in the 250 m case crosses the same Kelvin-Helmholtz pocket of unburned 
fuel (see Fig. [9]). Unlike in the 250 m case, however, oxygen fails to ignite (that is, 2 —f-+ 1? 
see Fig. [TTT) . and the high pressure spot that leads to a detonation in the lower resolution 
case eventually dissipates without initiating a detonation (see Fig. [10]). Thus in this case 
(which is the highest resolution simulation we were able to carry out), the compression wave 
leads to increased carbon burning in the region, but does not trigger an outburst of oxygen 
burning and is not sufficient to cause initiation of a detonation. This may be due to the 
fact that the temperature gradients, and therefore the induction time gradients, are much 
steeper, even in the cusp, compared to the lower resolution case (see fig. [8] bottom). 

Unfortunately, the computational demand of this simulation was great enough that we 
were unable to continue it to later times. At the time we stopped the simulation, the burning 
front at the head of the jet is still moving toward higher densities, and since the critical length 
scale for the initiation of a de t onati on is a strongly decreasing function of the density (see 



table 4 from ISeitenzahl et al.l (120091 )). a detonation may occur ahead of the jet at a later 
time. In addition, the strong sound wave that is involved in triggering the detonation in the 
500 m and 250 m cases is still moving toward the symmetry axis, which will increase the 
compression in the wave as it approaches the axis, making a detonation more likely. 



4. Discussion 



4.1. Comparison with Earlier Work 



Plewal ( 120071 ) states that in his GCD models the detonation a lways forms via a shock 



to detonation transition (SDT) (IBdzil fc Kapilalll992l ; ISharpell2002l ). In our results, the ini- 
tiation of the detonation is not tied to a shock. The 500 m and 250 m cases are somewhat 
similar in spirit to an SDT, since they also rely on the preconditioning of fuel by a hydro- 
dynamic wave. However, the wave is not a shock and the flows leading up to detonation are 
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not supersonic. While both simulations were performed using flash there are significant 
differences in the employed flame models, energetics, and initial conditions, which co uld ac- 
count for the disagreement. Another difference is the 8 km resolution employed in the iPlewa 
( 120071 ) study, which is 2-64 times coarser than in our study. 



4.2. Initiation of a Detonation in the GCD Model 

Initiation of a detonation occurs in three of the fours simulations we carried out. In 
all three cases, detonation occurs via the gradient mechanism. We consider the detonation 
in the 4 km case to be an artifact of the low resolution and numerical diffusion, which lead 
to a much shallower temperature gradient, and hence larger pre-heated region, than in the 
higher resolution cases. In the other two cases in which a detonation occurs, it is due to 
the passage of a compression front (500 m) or the passage of a compression front followed 
by a reflected wave (250 m) through the temperature gradient ahead of the jet, broadening 
the temperature gradient and making it shallow enough (and therefore the induction time 
long enough) that the gradient mechanism can operate. These mark the first simulations of 
gradient- induced detonations in a whole star simulation. 

In the fourth case (125 m), the highest resolution simulation we carried out, initiation 
of a detonation had not yet occurred by the time we stopped the simulation because of the 
computational demand it required. Whether a detonation eventually occurs in this case is 
an open question. We are therefore unable to state that initiation of a detonation via the 
gradient mechanism occurs in all of the higher resolution simulations. 

Our studies confirm that the temperatures and densities needed for initiation of a 
detonation to occur are readily obtained. The primary challenge is the necessity of pre- 
conditioning the fuel so that a shallow enough temperature gradient (and therefore a large 
enough induction time) occurs. Sufficiently shallow temperature gradients exist in the jet, 
but the material within the jet is mostly ash, due to the highly reactive nature of the fuel. 
Consequently, it is necessary that a sufficiently shallow temperature gradient occur in the 
unburned fuel ahead of the jet. For the resolutions explored in this paper, the spatial extent 
of the temperature gradient decreases as the resolution becomes finer, and the temperature 
gradient steepens with increasing resolution. It is therefore not too surprising that the 125 
m case (the highest resolution simulation we carried out) fails to detonate at the location 
where a detonation occurs in the 500 m and 250 m cases, or anywhere else (at least until 
the time the evolution was stopped). This is somewhat disconcerting. 



However, apart from the fact that an admixture of helium (e.g. lYoon fc Langerl 120051 ) 
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would lead more robustly to a detonation (see Table 12 in ISeitenzahl et al.l (120091 )). there 
are several reasons that suggest a detonation might still occur in the highest resolution case. 
First, the hot jet is still moving toward higher densities, and since the critical length scale 
fo r the initiation o f a det onation is a strongly decreasing function of density (see Table 4 



m iSeitenzahl et al.l (120091 )). a detonation may still occur ahead of the jet at a later time. 
Second, the strong sound wave that is involved in triggering the detonation in the 500 m and 
250 m cases is still moving toward the symmetry axis, which will increase the compression 
in the wave as it approaches the axis, making a detonation more likely. Third, the higher 
the resolution, the more turbulent shear flow structures emerge on small scales at the jet- 
fuel boundary (see fig. [12]). Before the high resolution truncated cone was imposed at 
t = 2.2124 s, the resolution of the simulation was 4 km, a scale at which Kelvin-Helmholtz 
shear instabilities are effectively suppressed. After the resolution in the truncated cone is 
suddenly increased, there is no structure initially on scales < 4 km. The turbulence on 
smaller scales requires time to develop; if given time to fully mature, it will smear the fuel- 
ash interface, leading to shallower induction time gradients, which will facilitate initiation 
of a detonation. 



In this regard, we note that iKnystautas et al.l (119791 ) and iThomas fc Jones! (120001 ) list 
three criteria that need to be fulfilled for a transition to a detonation in the turbulent shear 
layer of a jet: (1) large-scale energetic eddies in the unburned gas (2) strong fine turbulence 
sufficient to efficiently mix hot ash into the fuel; and (3) the establishment of an induction 
time gradient in the mixed region such that the SWACER mechanism can act. All three of 
these criteria appear to be fulfilled in our case. 



5. Conclusion 



The high-resolution simulations reported in this paper suggest that the detonations in 
lower resolution simulations (i.e., the 4 km case in this paper and the lower resolution simu- 
lations in earlier studies) are largely numerical. Two high-resolution simulations (500 m and 
250 m) demonstrate that pre-conditioning of a fuel pocket surrounded by Kelvin-Helmholtz 
billows by a compression wave moving through the pocket can lead to a gradient initiated 
detonation. For the resolutions explored in this paper, the spatial extent of the tempera- 
ture gradient decreases as the resolution becomes finer, causing the temperature gradient 
to steepen with increasing resolution. Consequently, the temperature gradient becomes too 
steep (and therefore the induction time becomes too short) in the highest resolution simu- 
lation (125 m) for a detonation to occur at the time and location at which it occurs in the 
500 m and 250 m simulations, and did not occur by the time we had to stop the simulation 
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because of its computational demand. However, it may detonate later. In addition, shear- 
induced turbulence, which surrounds the jet and is largely suppressed at low resolution, may 
broaden the temperature gradient, leading to conditions favorable for the initiation of a det- 
onation, especially if the influence of compression waves on the turbulent layer is considered. 
This suggests that future studies should be done (ideally in 3D) at even higher resolution, 
beginning earlier and running for longer times, in order to capture the development of the 
turbulence surrounding the jet more self consistently. 
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Fig. 1. — 500 m resolution. White line is p = 7 x 10 6 g cm -3 iso-density contour. X- velocity 
of material showing matter converging at the South Pole at time t = 2.2124 s, 54 ms before 
detonation. The accretion shock is clearly visible near x = 400 km off axis in this plot. 




Fig. 2. — 500 m resolution. White line is p = 7 x 10 6 g cm -3 iso-density contour. Y-velocity 
of material showing jet penetrating inward 54 ms before detonation. The internal shock see 
near y = —2200 km is not the location of the initiation. 
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Fig. 3. — 500 m resolution. White line is p = 7 x 10 6 g cm -3 iso-density contour. In 
this pressure plot of the collision region, 54 ms before detonation, the anatomy of the jet, 
including the various shocks and hydrodynamic waves, is clearly visible. 
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Fig. 4. — 500 m resolution. Left: Reaction progress variable 0i. Right: R eaction progress 



varia ble 02, a tracer for the progress of the burning to Si-group elements (ITownsley et al. 



20071 ) . Barely visible (near the center of the plots) is the detonation emerging at the boundary 
separating the carbon depleted jet region and the unprocessed fuel. 
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Fig. 5. — 500 m resolution. Blocks of fine resolution (16 x 16 zones per block, 500 m per 
zone) are used to cover a conical region near the tip of the inward moving jet, here visible 
as the region of high temperature resulting from the shock heating and subsequent carbon 
burning, 54 ms before detonation. 
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Fig. 6. — 4.0 km resolution. Left column shows a time sequence of the temperature. Right 
column shows rea ction progress variabl e <p2, a tracer for the progress of the burning to 
Si-group elements (ITownsley et al.l 120071 ). 
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Fig. 7. — 500 m resolution. Left: Detonation forms in compressed material. Right: Tem- 
perature in bulk of compressed region is not high enough for direct initiation of detonation. 



- 23 - 



time - 2,26654a 



time - 2.3S654s 



E 
u 




iuss a.1911 

i cm] 



time = 2,3«B7» 




O.KH 0.20C 



i [ltf cm] 



time = 2.26427s 




E 
u 



E 

41 



41U 0.190 
li [ICf cm] 



time = 2.a<jJB7i 



S.2Q6 DlKH 
■ [ltf cm] 



time = 2.26427b 




■3.0 a 




■Hi 



I [10" em] 



HZO* 




■4.05 



CZW 0.3 
» [ID 1 em] 



Fig. 8. — Temperature (left) and corresponding logarithm of the constant volume induction 
time Tin, (right) for 500 m, 250 m, and 125 m resolution (top to bottom). The sidelength 
of each plot is 6 km. r ijCV (see eq. [Q was set to 1 s for X$ > 0.9. For 500 m and 250 m the 
time of the snapshots is just before the detonation shock forms. For 125 m the time of the 
snapshots is just after the compression wave has passed the feature near x = 208 km. 
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Fig. 9. — Shown in each column is a time sequence of the density. Left: 250 m resolution. 
Right: 125 m resolution. 
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Fig. 10. — Shown in each column is a time sequence of the pressure. Left: 250 m resolution. 
Right: 125 m resolution. 
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Fig. 11. — Shown in each column is a time sequence of th e reaction progress va riable 02, a 
tracer for the progress of the burning to Si-group elements (ITownsley et al.ll2007l ). Left: 250 
m resolution. Right: 125 m resolution. 



-27- 



ilma = 2 time = P_.2652Sk 




0.1 0.2 5,5 0,4 0,1 O.Z G.3 0,* 

* [10 1 cm] K [1Q 1 cm] 



Fig. 12. — Comparison of density in the jet for 500 m resolution (left) and 125 m resolution 
(right). While the 125 m case shows turbulent structure on smaller scales, there is overall 
qualitative agreement. 
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Table 1. Detonation conditions 



Resolution 


Time 


Location (x,y) 


Fuel Density 


T 

1 max 


To 


[10 2 cm] 


[s] 


[10 8 cm] 


[10 7 g cm" 3 ] 


[10 9 K] 


[10 7 K] 


4,000 


2.204 


(0.00,-2.04) 


1.1 


-3.0 


-10 


500 


2.266 


(0.19,-1.95) 


1.2—1.6 


-3.0 


-20 


250 


2.264 


(0.21,-1.97) 


1.2—1.6 


-3.0 


-20 


125 


N/A 


N/A 


N/A 


N/A 


N/A 



